Spin-wave contributions to nuclear magnetic relaxation in magnetic metals 
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The longitudinal and transverse nuclear magnetic relaxation rates 1 /Ti (T) and 1 /T2 (T) are cal- 
culated for three- and two-dimensional (3D and 2D) metallic ferro- and antiferromagnets (FM and 
AFM) with localized magnetic moments in the spin-wave temperature region. The contribution 
of the one-magnon decay processes is strongly enhanced in comparison with the standard T-linear 
Korringa term, especially for the FM case. For the 3D AFM case this contribution diverges logarith- 
mically, the divergence being cut at the magnon gap ujo due to magnetic anisotropy, and for the 2D 
(N 1 AFM case as lu 1 . The electron-magnon scattering processes yield T 2 Iii(T/ujo) and T 2 /<x>q^ 2 -terms 

in 1/Ti for the 3D AFM and 2D FM cases, respectively. The two- magnon ("Raman") contributions 
are investigated and demonstrated to be large in the 2D FM case. Peculiarities of the isotropic 2D 
limit (where the correlation length is very large) are analyzed. 



> 

00 " 



e 



- 1—^ 

X 



a- 



I. INTRODUCTION 



Nuclear magnetic resonance (NMR), which is one of most powerful tools for investigating various physical properties, 
has a number of peculiarities for magnetically ordered materials. Last time, a number of new classes of magnets 
have been studied by this method, e.g., heavy- fermion compounds [E[, ferromagnetic films and monolayers low- 
dimensional systems including copper-oxide perovskites ||, etc. Thus the problem of theoretical description of various 
NMR characteristics of magnets is topical again. This problem was already a subject of great interest since 50-60s 
when the interaction of nuclear magnetic moments with spin waves in localized-spin Heisenberg model was studied 
^J,^|. However, this model is inadequate to describe the most interesting systems mentioned above where the role 
of conduction electrons is essential in magnetic properties. Usually the data on the longitudinal nuclear magnetic 
relaxation rate 1/Ti (this NMR characteristic is probably most convenient to compare experimental results with 
theoretical predictions) are discussed within itinerant-electron models such as Hubbard model or phenomenological 
spin-fluctuation theories. Ueda and Moriya |6|Q] calculated the dependences 1/T\(T) for weak itinerant magnets, 
' especial attention being paid to the paramagnetic region, and obtained strong temperature effects. Later this approach 
was extended to the two-dimensional case and extensively developed in connection with high-T c superconductors and 
related compounds ^|J9). On the other hand, in a number of systems (e.g., in most rare-earth compounds which are 
ON \ also a subject of NMR investigations, see, e.g., Refs. |h]]) the s — d(f) exchange model with well-separated localized 
ON ■ and itinerant magnetic subsystems is more adequate. Magnetic properties in such a situation differ essentially from 



those in the paramagnon regime (see, e.g., discussion in Refs. [g_2jl3 ). At the same time, the contributions to nuclear 
magnetic relaxation rate owing to electron-magnon interaction are not investigated in detail. 

In the present work we obtain the dependences of l/Ti(T) and the linewidth 1/T2(T) in the spin- wave region for 
three- and two-dimensional (3D and 2D) metallic magnets with well-defined local magnetic moments. In Sect. 2 we 
discuss the general formalism and physical picture of hyperfine interactions. In Sects. 3 and 4 we calculate various 
contributions to the relaxation rates in metallic ferro- and antiferromagnets. In Sect. 5 we analyze the isotropic 2D 
O ■ case where at finite temperatures the long-range order is absent, but the correlation length is very large. 



II. HYPERFINE INTERACTIONS 

We start from the standard Hamiltonian of the hyperfine interaction |l4| 

H hf = hi, h a = A al3 Sp (1) 

A being the hyperfine interaction matrix, which contains the Fermi (contact) and dipole-dipole contributions, 

A a(j = A F S a p + A d $. (2) 

The Fermi hyperfine interaction is proportional to the electron density at the nucleus and therefore only s-states 
participate in it, the contribution of core s-states (which are polarized due to local magnetic moments) being much 
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larger than of conduction electrons. It is just the consequence of considerably smaller localization area (and therefore 
higher density on nuclei) for the core states. 

The dipole contribution to Hhf can be represented as [Q 

H dip = «| [ i (/+5 - + rs +j _ l T z S ^ F (o) + 7+5+^(2) + 2(/-5+ + l+S z )F^} + h.c. (3) 
J 2 6 3 

where 

F (0) = ((1 -3cos 2 6»)/r 3 ),F (1) = (sin6>cos(9exp(-^)/r 3 ), 

= (sin 2 0cxp(-2^)/r 3 ),a = -\ 1&ln (4) 

Here (...) is the average over the electron subsystem states, j e and j n are gyromagnetic ratios for electron and nuclear 
moments, respectively. In the case of the local cubic symmetry we have F^ = 0. It is obvious that magnetic /- or 
d-electrons dominate also in dipole interactions because of large spin polarization. Hence the direct interaction of 
nuclear spins with that of conduction electrons can be neglected in magnets with well-defined local magnetic moments. 
Nevertheless, conduction electrons do effect nuclear relaxation via their influence on the local-moment system; besides 
that, as we shall see below, such contributions possess large exchange enhancement factors. The investigation of these 
effects is one of the main aims of this work. 

A general way to consider all these contributions is using the Green's function method which leads to the following 
expression for the longitudinal nuclear magnetic relaxation rate Jig] 



JY = -£ Im E<W = q>>">>- ( 5 ) 

q 



(w„ = (h z ) < T is the NMR frequency). As follows from 



Then we derive 



h- = (A F + ioF<°>)£- + aF (2 \S+ + 2aF^S z , (7) 

h* = (A F - -aF ( V)S z + a(F^S + + aF^*S~) (8) 
3 



^ = l {[(AF + r F(0))2+fl2|F(2)|2] ^ 

+2a(A F + -aF^)ReF^K++ + Aa 2 \F^\ 2 K zz } (9) 
3 



— = — + -{(A F - -aF^fK" + a 2 [2\F^\ 2 K+- + (F^) 2 K++ + (F^*) 2 K ]} (10) 
T2 2Ti 2 3 

where the quantities K a/3 are defined by 

K aP = -- lim Im^((^|^ q )) w / W (11) 
q 

As we shall see below in Sect. 5, ui n ^ which enters (^|) (but not the second term of ([lO])) may become important in 
the case of very small magnetic anisotropy. 

Formula (|5|) has a rather general character. On the other hand, the problem of calculating the NMR linewidth is 
much more complicated. The Moriya formula (^|) is in fact applicable only in the case where the line has the Lorentz 
form (i.e., charactersitic frequency of hyperfine field fluctuations is large in comparison with their amplitude) [ pT[ . 
In insulating crystals the latter condition is usually violated, the lineform being close to Gaussian with the width 
determined by dipole interactions of nuclear spins. At the same time, in metals the Korringa relaxation described by 
the formula (||) usually dominates, so that we will use this. A peculiar case is provided by conducting systems which 
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are on the borderline of ferro- or antiferromagnetic instability (i.e. with large correlation length £), e.g., copper-oxide 
superconductors [fll| . Under this condition, the anisotropic Ruderman-Kittel interaction between nuclear spins turns 
out to be greatly ehanced and dominates over the dipole interaction. The lineform turns out to be Gaussian with the 
width being estimated as 

(^) 2 «A 2 £ x 2 ( q ,c = o) ( 12 ) 

^ 2 ' q 

where x(q, ^-0 is the dynamical spin susceptibility of electron system. We will use this result in Sect. 5 when discussing 
2D systems which do possess large correlation length. 



III. FERROMAGNETIC METALS 



We proceed with the s — d(f) exchange model Hamiltonian 

H = ^ *kCkff C k<7 - I ^2 Si<T a pc\a C iP + ^2 ^q S -q S q + H a (13) 
ko- iafi q 

where is the band energy, Si and S q are spin-density operators and their Fourier transforms, a are the Pauli 
matrices, H a is the anisotropy Hamiltonian which results in occurrence of the gap loq in the spin-wave spectrum. For 
convenience we include explicitly in the Hamiltonian the Heisenberg exchange interaction with the parameters J q , 
although really this may be, e.g., the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction. It should be noted that 
similar results may be reproduced for the localized-moment Hubbard magnets (cf. fl^ , |l6|| ). 

First we consider the ferromagnetic (FM) case. Then K ++ = and the relaxation rates (^), JlO| ) are the sums of 
transverse (oc ) and longitudinal (oc K zz ) terms. Passing to the magnon representation we obtain 

((^|SZ q )) w = 2S/[w - c q + i 7q (w)] (14) 
where w q = 2S( J q — Jo) + loq is the magnon frequency, 7 q (w) oc w is the magnon damping. Then we have 

K+- =2SY^4 (15) 

(cf. Refs. (T^JlTl l). The damping in the denominator of (|T~5 |) can be neglected for both localized-moment and itinerant- 
electron magnets (in the latter case the expression (|14|) corresponds to the RPA structure, see Ref. fl2|| ) due to 
smallness of uj n . On the contrary, temperature dependences of magnetization, resistivity etc. in weak itinerant 
magnets are just determined by the damping in the denominator, i.e. by paramagnon excitations rather than by spin 
waves fjj. 

The damping owing to the one-magnon decay processes is given by the well-known expression 

7 £V) = -2^T 2 S]T(n kT - n k _ qi ) (16) 

k 

X6(UJ + tk| — ^k-ql) — 27r/ 2 StoXq 

where tk<r = ^k — crIS, is referred to the Fermi level, = n{t-^ a ) is the Fermi function, 

Aq = $>(*kTW*k-qi)- (17) 

k 

The linearity of spin fluctuation damping in u> is the characteristic property of metals. According to (^) this leads 
to T- linear contributions to 1 fT\ which is the Korringa law (Ti| . Note that the simplest expression for the Korringa 
relaxation 

1/Ta ~ 1/T 2 ~ A 2 PWi T, (18) 

where A is an effective hyperfine interaction constant, p a are the partial densities of electron states at the Fermi level, 
is practically never applicable for magnetic metals (e.g., exchange enhancement factors can change even the order of 



3 



magnitude of 1/Ti [[Tjjl 7|] ) . Accurate expression for the "Korringa" contribution in the case under consideration can 
be derived by the substitution ( |l5|) and (|l6| ) into (||). 

The damping ( |l6| ) has the threshold value of q, which is determined by the spin splitting A = 2|/|5, q* — A/vf 
(vf is the electron velocity at the Fermi level). The quantity q* determines a characteristic temperature and energy 
scale 

OJ* = uj{q*)=V{q*) 2 ~(A/v F ) 2 T c (19) 

with V the spin-wave stiffness. 

Besides 3D magnets, consideration of the 2D case is of interest (this may be relevant, e.g., for layered magnets and 
ferromagnetic films; for more details see Sect. 5). We have 

A q = %- 9 >{^ AV / 2 j^ (20) 

After integration for the parabolic electron spectrum (q* plays the role of the lower cutoff), the one-magnon damping 
contribution to (|l5|) takes the form 

^) K+ - = mL x Nf £ = 3 (21) 

V 2 m 2 \l/(irq*), D = 2 y ' 

with 

_mfi„ (k Fa /iT, D = 3 

Pa -^T X \1, D = 2 [22) 

m the electron effective mass, Slo the lattice cell volume (area). Thus in the 3D case the factor of I 2 is canceled, and 
the factor of I~ l occurs in the 2D case and we obtain a strongly enhanced T- linear Korringa- type term (remember 
that V ~ J ~ I 2 p for the RKKY interaction). This means that the contribution of conduction electrons to T- 
linear relaxation rate via their interaction with localized spins is indeed much more important than the "direct" 
contribution (|l8|): perturbation theory in the s — d exchange coupling parameter / turns out to be singular. Earlier 
such contributions (for the 3D case) were calculated by Weger [Q and Moriya ^] for iron-group metals. However, 
Moriya has concluded that for these materials they are not important in comparison with orbital current contributions. 
In the case under consideration (where magnetic subsystem is well separated from the conductivity electrons) the 
situation is different and the spin- wave contribution in 1/Ti is normally the most important. 

The one-magnon decay contribution ( pl| ) is absent for so-called half- metallic ferromagnets, e.g., some Heusler alloys, 
where electron states with one spin projection only are presented at the Fermi surface pl|pT[ . In such a situation we 
have to consider two-magnon scattering processes. In this connection, it is worthwhile to note an important difference 
between relaxation processes via phonons and via magnons. The main difference is due to the gap in magnon 
spectrum. Usually ojq > oj n and therefore one-magnon processes contribute to the relaxation rate due to magnon 
damping only (cf. discussion of the phonon-induced relaxation processes in Ref. [flif ). However, the mechanisms of 
magnon damping in magnetic dielectrics (magnon-magnon interactions) are different from those in magnetic metals 
and degenerate semiconductors |p2}[l3| . 

The damping in a conducting ferromagnet owing to electron-magnon (two-magnon) scattering processes is calculated 
in Refs. f^fjlill and has the form 

S 2 )i, ,\ / + . +. \ 2 Vp 



f — ' V tk+a — + 2a IS J <9£k dcu u 

kper ^ ' f 



x5(£ k -<k-p+q) (23) 
where iV p = N(uj-p) is the Bose function. Substituting this into ( |l5|) and performing integration we obtain for D = 3 



128TT 2 Sm 2 W 2 \8M 3 u;* : T » lo* 

where C,(z) is the Riemann function, 



Ms = / dx 



1 x exp x 



,.2 



(exp x 2 — l) 2 



0.65 (25) 
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The contribution ( p4| ) should play the dominant role in the half-metallic ferromagnets [jl7). Besides that, this contri- 
bution may modify considerably the temperature dependence of 1/Ti in "usual" ferromagnets, a crossover from T 5 / 2 
to T 3 / 2 dependence of the correction taking place. 

1 /2 \ K 

For D — 2 at T, lo* ^> loq small magnon momenta of order of (ujq/T>) 1 make the main contribution to (15). To 
calculate the integral one can use the high-temperature expression for N p — T/aj p . As a result, one gets 

6(2)R+ _ = nlk F M 2 



with 



M 2 



tt/2 

dx f dip sin 2 <p 



o 1 + x 2 J ( sm 2 ^ + x 2\ 3 / 2 







dy 



1.23 (27) 



Thus in the 2D FM case, in contrast with 3D one, the relaxation rate 1/Ti is strongly dependent on the anisotropy 
gap. 

Consider now the second term in the tranverse relaxation rate l/Ta(T) (|l^), which is normally determined by 
K zz , and the longitudinal contribution to relaxation rate 1/Ti in @, which is due to dipole-dipole interactions with 
the characteristic constant A ~ a\F^\. The simplest calculation from the longitudinal Green's function for the 
localized-spin subsystem gives 

W'q»"=E , (28) 

qp ^ p ' 

The quantity ( p9| ) has been considered in Refs. p3|,p[ as a contribution to the NMR line width 1/T 2 . The integration 
in the 3D case gives the logarithmic singularity 

K zz = ° Tin — (30) 

16tT 4 T> 3 luo V ' 

For D = 2 this singular term is inversely proportional to the magnetic anisotropy parameter and very large: 

For small enough uiq and A ~ yl this contribution can dominate over the "Korringa" contribution ( pll) in 1/Ti at 
T > u>o/\Ip\. The leading contribution to K zz from the s — d interaction is determined by 



<5«S q |^ q )) w = 2/ 2 S]T 



x nkjXj- ~ n k + p - qT ) + N p (n kl - n k +p-qt) / 32 s 
tkl — ik+p-qT + crw — o; p 

However, it is not singular in ujq and practically never important. 

IV. ANTIFERROMAGNETIC METALS 

Now we consider the spiral antiferromagnetic (AFM) structure along the x-axis with the wavevector Q 

(5?) = S*cosQR,, (Sf) = S*sinQR 4 , (Sf) = 
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We introduce the local coordinate system 

S? = Sf cos QR, - Sf sin QR 4 , 

Sf = Sf cos QR 4 + S? sin QR 4 , Sf = Sf 

Further we pass from spin operators Sj to the spin deviation operators b\,bi and, by the canonical transformation 
= u q/^q — w q/5-qj to the magnon operators. Hereafter we consider for simplicity two-sublattice AFM ordering (2Q 

is equal to a reciprocal lattice vector, so that cos 2 QR, = 1, sin 2 QR, = 0). 

Calculating the Green's functions to second order in I (to second order in the formal quasiclassical parameter 1/2S, 

cf. Refs. fl3|]24|| ) we derive 

(M^», = 7 r + (33) 

W-M)" = 7 r v °2r VTJvT (34) 

with 

Cquj = S(Jq^ q |ta) + J q °* — 2Jqq) + y^[Cp$pq^ 

P 

-(C p - D p )<P pQ0 + </>+ qw + </> pq J + Sq (35) 

^ qu = = S(J qw — ^Q-(- q , w ) + y ] Dp^pquj + hq 

P 

The s — d exchange contributions of the first order in 1/2S correspond to the RKKY approximation 

_ - Hk-q 

ii 

k 



(nk = ^(ik) is the Fermi function), the second summand in (pq ) being the w-dependent RKKY indirect exchange 
interaction. The function <i>, which determines the second-order corrections, is given by 

$pqa, = (</>p qa) - ^pqj/^p, (37) 
, ± _ T 2 "k(l - rtk+p-q) + iVXjgpKnk - »k+p-q) 



pq^ 



= i 2 Y 



where 

uj p = (C 2 p - Dl) 1 ' 2 = [4S 2 (J P - J Q )(J Q+P - Jq) +c 2 ] 1 /2 

is the magnon frequency to zeroth order in /and 1/2S. The ^-independent corrections <? q , fo q that describe the 
"direct" magnon-magnon interaction are written down in Refs. [ , p4j , 

Now we consider the effects of electron-magnon interaction. The intrasubband one-magnon damping (which is 
absent in the FM case) is finite at arbitrarily small q J25|. Similar to the FM case, the contributions of intersubband 
transitions (which correspond to small |q — Q|) are cut at the characteristic temperature and energy scale 

uj* = uj(q*) = cq* ~ (A/v F )T N , (38) 

We have 

2S 

K + - = lim ImVw-'Cau^, (39) 

q 

and the term with K ++ in (||) vanishes due to the property -D qw = — -D q +Q W . The one-magnon contribution owing to 
the imaginary part of ( |36|) in the 3D case takes after integration the form 

^-^f^ln^+P^^V (40) 
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Here c is the magnon velocity defined by Wp = ^p + Q = + c 2 p 2 , 

P p = I 2 lim |q - p| V 5(t k )<5(i k _ q+p ), (41) 

q^O * — * 
k 

(the quantity Pq depends, generally speaking, on the direction of the vector q, see Refs. |?],^6],^7|] ) , the second logarithm 
in the brackets of ( fhi| ) contains the cutoff 

uj* = ^ 2 + {u*y 

The "enhancement" factor in (^0|) is smaller than in the FM case because of the linear dispersion law of magnons, 
but this contribution still dominates over the "usual" Korringa term (|l8|). Besides that, a large logarithmic factor 
occurs (in the isotropic case, this is cut at u n only). Note that a similar logarithmic singularity in 1/Xi takes place 
for 3D itinerant-electron antiferromagnets ^] . It is interesting that the intersubband contribution does not lead here 
to enhancing the singularity, unlike the situation for the magnon damping, magnetic and transport properties p8| , |l3| . 
Under the "nesting" conditions (tk+Q — — ^k in a large part of the Fermi surface) the singularity is not enhanced as 
well. 

The singularity becomes stronger in the 2D case where integration gives 

S W K+ - = (l Po + Pq arctan ^) (42 ) 

TTCLUo V 2 UJ* / 

This fact may be important when treating experimental data on layered AFM metals. 

The contribution owing to electron-magnon scattering processes is determined by the imaginary part of the function 
([[t]). After a little manipulation we obtain 

S^K+- ~ 2SL V -4- (-^) [Po + Pq?(?)] (43) 



where L — 2S(Jq — Jq), 4>{q < q*) = 0, cf>(q q*) = 1. The integration in the 3D case yields 



s(2)R+ _ = [PomuJo) + P Q /(T, wS)] (44) 



where 



duj 



wo \ w 2 uj 
Thus we have 1 jT\ oc T 2 In T. In the 2D case we derive 



Tin— fln^-iln— |,T»u,„ (45) 



so that the singularity is not enhanced in comparison with the 3D case. 

The contributions owing to longitudinal fluctuations will be estimated for the localized subsystem only. We obtain 

L2 ( dN P 



-af ]*K-p) (47) 



The corresponding contribution to I/T2 was considered in Ref. [Q. The term in the longitudinal relaxation rate 
determined by (^) is estimated as 

Provided that the dipole-dipolc contributions in @ are considerable (^4 ~ A), this term can dominate over the 
"Korringa" term @ of order of A 2 I 2 p 2 T\n \ J/uj \/J 2 at T/\J\ > \Ip\\n 1/2 \ J/u \ only. Note that this two-magnon 
contribution is similar to the two-phonon (Raman) contribution in the spin-lattice relaxation. The existence of the 
gap Cl>o is not important here (at least if it is sufficiently small) , but the matrix elements of interaction of nuclear spins 
with magnons are singular, unlike those for acoustic phonons (|M q ^o| 2 ~ l/<7 instead of q ). Therefore we have a T 3 
law instead of T 7 one for the phonon scattering [|l4| . 
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V. ISOTROPIC 2D CASE AND NMR IN LAYERED AND FRUSTRATED MAGNETS 



Now we investigate the case of layered magnets, in particular, the isotropic 2D limit. A detailed treatment of the 
spin correlation functions and corresponding spin-fluctuation contributions to 1/T± in the isotropic 2D Heisenberg 
antiferromagnets with w„ —> was performed in Ref. p9] j . Here we calculate also corrections owing to electron- magnon 
interaction. 

The magnetic ordering temperature is determined by magnetic anisotropy or interlayer coupling, 

T M ~ |J|S 2 /M|J|S7max{^ , |J'|}) (49) 

(for more details see, e.g., Refs. |50|). Despite the absence of the long-range ordering (LRO) at finite temperatures, 
spin-wave description holds even in the pure 2D isotropic case in the broad temperature region up to T ~ | J\S (i.e., 
Tm \ J\S 2 ) owing to strong short-range order (SRO). In a more general case of finite ujq and J', this description holds 
at T ^ Tm- A possibilty to describe LRO without introducing anomalous averages (like sublattice magnetization) in 
terms of singularities of the spin correlation function was demonstrated in Ref. [pl| . Such an approach enables one to 
obtain an unified description of ordered and disordered phases. In the pure 2D case the magnetization (or sublattice 
magnetization for the AFM case) S is replaced in both magnetic and electronic properties by the square root of the 
Ornstein-Cernike peak (see, e.g., p3]). The gap in the effective spin-wave spectrum appears at finite temperatures, 
which is determined by the inverse correlation length. The correlation length in the situation under consideration is 



estimated as 1 33 



£ oc exp (tt| J\S 2 /T) (50) 

As shows the two-loop scaling theory f34j| , the preexponentional factor is temperature independent; quantum effects 
can renormalize the exchange parameter J fl3Cfl . 

To describe formally NMR in the absence of LRO ((Si) — 0) we follow to Ref. |H| and consider the autocorrelation 
function of the nuclear spin I ]3q |. Performing calculations with the simplest Hamiltonian H^f = AISi to second 
order in A we derive 

(7+,7-r = ^+i)[-w+s(w)] 

with the memory function 

f°° 1 
E(w) = A 2 / dteM^t) X;<5* q (*)5* + ^SZ^{t)S+), (51) 

q 

r°° 1 

(S^(t)S?)= / du J eM^t)jf(tu),jf(u J ) = --N(u J )Im((SP\S« cl }} UJ . (52) 

J — oo f 

As discussed in Ref. the spectral density J^iw) contains an almost singular contribution 

5jf(u) cx A q _ Q A w (53) 

where A q and A u are delta-like functions smeared at the scales q ~ and u> ~ u>^ with the characteristic spin- 
fluctuation energy cu^ ~ D£~ 2 (FM), ui£ ~ c^ 1 (AFM). To obtain the singular term in with the correct factor 
of S 2 we can introduce a very small magnetic anisotropy (which does not violate time- reversal symmetry), so that 
the whole singular contribution passes to K^ z (cu) (cf. Ref. Then the term iA 2 S 2 /u occurs at w > uj^, and we 
obtain the expression 

(^•^-^^"(^as + JTs) (54) 

which describes precession of the nuclear spin with both frequencies ±AS. We see that the resonance picture holds at 
uj n ^> lo^ only. In the opposite case the NMR line is smeared, but we can calculate the quantity 1/Ti according to 

Provided that ujq -c lu^ the quantity to^ plays a role of the gap in the magnon spectrum. Therefore at uo n <C 
the cutoffs in the singular contributions to 1/Ti described by (|42|), (E^), (0) are determined by very small inverse 
correlation length, so that they have very large values and possess unusual temperature behavior. 
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In the 2D FM case the expression ( prf ) is also applicable, but with another expression for loq^loo — > TJ>£ 2 which is 
exponentially small. We have 

<5 (2) (1/Ti) oc / 2 A 2 T 2 /^ 5/2 ^ 1/2 oc lW n HT 2 /V 3 . (55) 
In the isotropic 2£) AFM case we obtain from ( |42| ) 

<S (1) (1/Ti) oc 7 2 A 2 T/c^ oc I 2 u 2 £T/c 2 . (56) 

As follows from ©, @, 

<S (2) (1/Ti) oc 7 2 ^ 2 (T 2 /c 4 )ln 2 £ ~ const(T). (57) 

Possibility to observe such dependences experimentally is of great interest. Unfortunately, most experimental data for 
layered compunds deal with copper-oxide systems which are on the bordeline of AFM instability. The latter results 
in a specific temperature dependence of the spin susceptibility and strong deviations from the Korringa law. This 
makes separation of one- and two-magnon contributions impossible. 

At u> n 3> ^£ the cutoff in the singular contributions to 1/Ti described by (42), jig) , ( [jij ) is the NMR frequency u> n . 



However, such a cutoff is absent for the corresponding terms in I/T2 owing to the interaction A ~ al-F^). Provided 
that A ~ A we reproduce for these terms the dependences ([35j)-(^37|). However, as discussed in Sect. 2, for large £ the 
lincform turns out to be Gaussian owing to strong Ruderman-Kittel interaction between nuclear spins, the Moriya 
formula (^) is inapplicable, and one can estimate the linewidth from (|l2|). In the 2D case we obtain I/T2 oc £. A more 
detailed discussion of this situation and an application to copper-oxide systems is given in Ref. fTT| . 

We see that NMR investigations can be in principle used to obtain the temperature dependence of the correlation 
length. When crossing the magnetic ordering point in layered systems the NMR picture should not change radically. 
Formally, as follows from fl49|), ( |50| ) at T ~ Tm we have ln£ ~ In | JS 2 / max{wo, |</'|}|, so that the cutoffs are joined 
smoothly. 

A similar situation can take place for other systems with suppressed LRO and strong SRO, e.g., for frustrated 3D 
magnetic systems ||(| (where £ is also large, the lineform is Gaussian, and we obtain from ( |l2| ) I/T2 oc C 1 / 2 ). This 
may explain why the problem of detecting long-range magnetic ordering is frustrated systems with small ordered 
moments with the use of the NMR method so difficult. Indeed, the NMR data for heavy- fermion systems are doubtful 
and contradict to results of other experiments yj] . 



VI. CONCLUSIONS 



In the present paper we have investigated in detail various mechanism of nuclear magnetic relaxation in metallic 
ferro- and antiferromagnets in the spin- wave temperature region. In the most cases the main contribution to 1/Xi is of 
Korringa type, but its physical origin is more complicated than in paramagnetic metals. Formally, it results from the 
interaction of nuclear magnetic moments with the localized electronic subsystem with taking into account the "Stoner" 
(Landau) damping of spin waves via conduction electrons. This contribution is greatly enhanced in comparison with 
the standard Korringa term by inverse powers of exchange interaction (s — d(f) parameter), especially in ferromagnets. 
In 31? antiferromagnets such a contribution contains the logarithmic singularity which is cut at the gap in the magnon 
spectrum (magnetic anisotropy energy) lu . Thus we can conclude that the "Korringa" relaxation rate in magnetic 
metals should be much larger than in paramagnetic ones where the relaxation is determined by direct interaction 
with conduction electrons (such a term is also present in the magnetically ordered state, but is much smaller than 
the contribution discussed). In the 2D AFM case we have 1/Ti oc uJq 1 - In the isotropic limit the singularity in 1/T± 
is cut at very small inverse correlation length, so that the one-magnon contribution becomes very large and possesses 
unusual temperature behavior. 

Besides that, we have calculated contributions from more complicated magnon damping processes (electron-magnon 
scattering). For antiferromagnets and 2D ferromagnets they contain singular logarithmic or (in the 2D FM case) 
power-law factors which are also cut at ujq. These contributions may result in considerable deviations of the temper- 
ature dependence of 1/Ti from the linear Korringa law. In the 3D FM case this contribution is also noticeable and 
probably can be separated when fitting experimental data. For half-metallic ferromagnets, where the "Stoner" damp- 
ing is absent, this scattering should be the main nuclear relaxation mechanism (see the discussion of experimental 
data in the review p7|). 

Provided that the "longitudinal" matrix elements of dipole-dipole hyperfine interactions in (j^) are not too small, 
the two-magnon ("Raman") relaxation processes may be also important in 1/Ti, especially for 2D ferromagnets. 
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